Intermittency and universality in fully developed inviscid and weakly compressible 

turbulent flows 
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We performed high resolution numerical simulations of homogenous and isotropic compressible 
turbulence, with an average 3D Mach number close to 0.3. We study the statistical properties of 
intermittency for velocity, density and entropy. For the velocity field, which is the primary quantity 
that can be compared to the isotropic incompressible case, we find no statistical differences in its 
behavior in the inertial range due either to the slight compressibility or to the different dissipative 
mechanism. For the density field, we find evidence of "front-like" structures, although no shocks 
are produced by the simulation. 
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Fully-developed three-dimensional turbulence is char- 
acterized by an intermittent energy flux from large to 
small scales. For incompressible and viscous turbulent 
flows, dissipation occurs close to the Kolmogorov scale 
r\. According to the Kolmogorov theory, the statistical 
properties of turbulence are scale-invariant within the 
inertial range r\ <C r <C Lq, where Lq is the scale of 
energy forcing. Intermittency spoils dimensional scale 
invariance and is the origin of anomalous scaling 
There exist a number of numerical and experimental re- 
sults aimed at investigating the intermittency in turbu- 
lence for incompressible and viscous flows. On the other 
hand, few investigations have been reported so far for 
the case of weakly compressible and inviscid turbulence, 
relevant to many astrophysical and geophysical problems 
0, Q. In this letter we present a high-resolution nu- 
merical simulation of homogeneous and isotropic, three- 
dimensional compressible and inviscid turbulence. Our 
main purpose is to investigate intermittency in the iner- 
tial range and compare our finding against known results 
for incompressible viscous turbulence. For the case of 
driven and decaying supersonic turbulence similar prob- 
lems have been addressed in [1, For this simula- 
tion, we use the FLASH 3 component-based simulation 
framework. While the FLASH framework was primar- 
ily designed to treat compressible, reactive flows found 
in astrophysical environments [f|, it is generally appli- 
cable to many other types of fluid phenomena. For this 
simulation, only the compressible hydrodynamics module 
based on the higher-order Godunov Piecewise Parabolic 
Method (PPM) was used Q ■ The algorithmic methodol- 
ogy of the FLASH framework and the computer science 
aspects of this turbulence simulation have been described 



in further detail elsewhere [8(] . The equation of motions 
solved are the Euler fluid equations with forcing: 



dp 
di 



V ■ (vp) = 



dpv 

~dT 

dpE 



V • (vvp) = -VP + F 
V • [v(pE + P)\ = 



P = ( 7 - l) P U , E = U + -pv 2 



(1) 
(2) 
(3) 
(4) 



where p is the density, v the velocity, P the pressure, E 
the total energy, U the internal energy and 7 is the ratio 
of the specific heats in the system, equation (0} being the 
equation of state for our system. The effect of the large 
scale forcing F in ([2]) gives rise to a turbulent flow whose 
energy is transferred from scale Lo towards small scales 
Q. The energy input J vFd 3 x produces an increase of 
the internal energy U, which grows in time. One can 
easily show that the quantity P<9;V; represents the en- 
ergy transfer from kinetic to internal energy of the flow, 
which acts primarily on the smallest scales. The mean 
sound speed increases slightly in time as well, though 
the 3-D RMS Mach number is roughly 0.3 (1-D Mach 
number 0.17) throughout. The numerical simulation was 
done for isotropic and homogeneous forcing Q with a 
state-of-the-art resolution of 1856 3 grid points. The inte- 
gration in time was done for 3 eddy turnover times after 
an initial transient evolution beginning from rest. A key 
feature of the numerical simulation is that while the flow 
is formally inviscid, some viscosity is introduced as a re- 
sult of the numerical scheme. The numerical treatment 
of turbulent flows used in the FLASH simulation is some- 
times referred to as an Implicit Large Eddy Simulation 
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FIG. 1: Local slopes of the 8th-order (top) and 4th-order 
(bottom) longitudinal (blue lines) and transverse (red lines) 
structure functions as obtained from the present ILES (bold 
lines) and compared against a DNS of the Navier-Stokes equa- 
tions at comparable Reynolds numbers (narrow lines). 




FIG. 2: Two dimensional section of the density field p at a 
given time: large regions with smooth density variations are 
separated by sharp cliffs. 



(ILES), to be distinguished from a full Direct Numer- 
ical Simulation (DNS) of the Navier-Stokes equations. 
Thus, one might expect that the dynamics of the flow 
may significantly differ from incompressible and viscous 
high Reynolds number turbulence. 

A key question is therefore whether the statistical 
properties of the inertial range are significantly differ- 
ent from to the homogeneous and isotropic turbulence 
observed in Navier-Stokes. To answer this question, wc 
employ the concept of anomalous exponents in measur- 
ing the statistical properties of velocity fluctuations in 
the inertial range. In particular, we use the well known 
method of Extended Self-Similarity (ESS) which 
allows us to accurately estimate the anomalous expo- 
nents of the longitudinal = {[5v(r) ■ r] p ) (where 
5v(r) = v(x + r) — v(x) and (...) means average over 
the volume and in time) and transverse structure func- 
tions S [ P = {[Sv(r T )] p ) (where r T • v = 0). We de- 
note the corresponding scaling exponents by Cp{L) and 
C P (T). Our numerical result for ( p (L),(p(T) agrees re- 
markably well with previous numerical and experimental 
data. This is shown in figure [TJ which compares the ESS 
local slope d\og(Sp) / d\og(S2) , for both longitudinal and 
transverse structure functions p — 4 and p = 8 with the 
DNS simulations performed for incompressible Navier- 
Stokes equation at comparable Reynolds numbers [Io| . 
Figure Q] shows several interesting features. First, there 
exists a range of scales (r/r/ > 50) where the local slope is 
almost constant, i.e. where we can detect accurately an 
anomalous scaling exponent. Second, as one can see, in 
the inertial range our ILES results give exactly the same 
anomalous scaling as that obtained for the incompressible 
Navier-Stokes equation [lfj. There is a clear difference 
between our results and the Navier-Stokes case in the 



dissipation range (r/77 < 20), where the Navier-Stokes 
solutions show a well defined "dip" effect (as qualita- 
tively predicted by the multifractal theory [l2|), while 
ILES behaves differently. The different behavior in the 
dissipation range is expected since the ILES does not dis- 
sipate energy in the "standard" Navier-Stokes way. On 
the other hand, the remarkable agreement in the inertial 
range allows us to claim that the inertial range properties 
are independent of the dissipation mechanism. This is an 
important statement, which has been questioned several 
times in the past, and supports the conjecture that the 
statistical properties of turbulence in the inertial range 
are universal and independent of the dissipation mecha- 
nisms. This is one of our main results. 

We also note that the difference between longitudinal 
and transverse scaling exponents observed in homoge- 
neous and isotropic DNS [lfj, is also seen in our numer- 
ical results. This discrepancy is an open theoretical is- 
sue, not explainable using standard symmetry argument 
in homogeneous and isotropic turbulence [13|. Whether 
this remains true at higher Reynolds numbers is an open 
question (see also Q for a discussion on this point). 

Although the integration is formally inviscid, there is 
a net energy transfer from the turbulent kinetic energy 
l/2pv 2 to the internal energy. Thus we may consider 
that an effective viscosity v el{ is acting on the system. In 
order to estimate u B{{ we can proceed as if the Kolmogorov 
equation -with effective parameters- applies to our case: 

si L \r)=~e el[ r + 6 K ^Si L) (r) (5) 

A fit of our data with this formula gives, e eff = 
0.054, j/ cff = 8.3 • I0 -6 , which corresponds to a Kol- 
mogorov scale, 77 = (^ ff /eeff 4 )i equivalent to roughly half 
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FIG. 3: Compensated structure functions of the density 
D§(r)r~ z °° (solid line) and D$(r)r~ z °° (dotted line). In the 
insert, we show the scaling exponents of the density structure 
functions (•) and the entropy structure functions (a). Er- 
rors include both statistical fluctuations and the uncertainty 
in the fit to the inertial range. 



grid cell and to R\ ~ 600. The dynamical effects of 
the effective viscosity are however different from what 
one usually observes in the Navier-Stokes equations; i.e., 
the dissipation range does not behave the same as in the 
Navier-Stokes solutions. Thus, while we still observe an 
effective dissipation in the energy flux (i.e., in the third- 
order structure functions), the behavior in the dissipation 
range is different because the dissipation mechanism in 
our simulation is different: the increase of internal energy 
is due to viscosity resulting from the numerical scheme. 

In figure [2] we show the density field in the system at 
a given time. One can easily recognize the existence of 
large density gradients due to compressibility. Another 
interesting quantity to look at is the entropy S defined as 
S = \og(P/ p 1 ). Using equations ©-(HI) one can obtain 
the following equation for the entropy: d t S + v ■ VS = 0. 
This equation tells us that S satisfies an equation sim- 
ilar to the case of a passive scalar advected by the ve- 
locity vector v, although S cannot be considered a pas- 
sive quantity in this case. Strong variations are also de- 
tectable in the field S (not shown). Recently, the statis- 
tical properties of density fluctuations have been inves- 
tigated for supersonic turbulence characterized by large 
Mach number. A rather large effect is caused by the for- 
mation of shock waves and fronts 0, H, EH- Here, the 
3D Mach number is of order 0.3 on average and, conse- 
quently, we should expect that shock waves are not im- 
portant (although the compressibility degree may reach 
a maximum excursion where Mach ~ 0(1)). Actually, it 
has been shown in experiments and direct numerical sim- 
ulations of passive scalars, that front-like structures are 
Til . 17 1 - A front-like structure on a 




FIG. 4: Rescaled probability distribution function, 
P[5p(r)]r~ z °° . Different curves corresponds to different 
values of r/n in the inertial range (respectively, from top to 
bottom, 68,112,184, 288.) 



Q(x + r) — Q(x) ~ const for x + r and x selected on the 
two different sides of the front. If these "front-like" struc- 
tures play a significant role in the statistical fluctuations 
of Q, one should observe quite intermittent anomalous 
scaling for the structure functions of Q. In particular the 
anomalous exponents should approach a constant value 
for p — > oo. 

In order to study the statistical properties of the den- 
sity and entropy fluctuations, we can introduce the den- 
sity structure functions D p (r) = ([5p(r)] p ) and the en- 
tropy structure functions E p {r) = {[8S{r)f). The ILES 
result shows that Dp{r) and E p (r) are scaling func- 
tions of r in the inertial range; i.e., D p (r) ~ r Zp and 
E p (r) ~ r hp . The values of z p and h p are shown in 
the insert of figure [3] together with the compensated plot 
D$r~ Z5 and D§r~ Zti . From figure [3] one can appreciate 
the quality of the scaling in the inertial range. We re- 
mark that the scaling exponents z p and h p show quite 
anomalous behavior. For p < 3 the values of z p are larger 
than the corresponding values of ( P {L) and ( P (T) which 
means that density is somehow smoother than the veloc- 
ity field. The anomalous exponents for both the density 
and the entropy structure functions become constant at 
large order. In particular, defining the saturation expo- 
nents to be Zoo and hoo, we estimate Zoo = 1.62 ± 0.10 
and hoo = 1-00 ± 0.10. A similar analysis for the pres- 
sure field (not shown) shows that the structure func- 
tions of the pressure field P have the same scaling ex- 
ponents and of the same saturation exponents as those 
for the density field. Using equation (Q} and assuming 
stationarity, one can derive an exact relation for the "en- 
ergy" flux of density fluctuations ((Sp(r)) 2 Sv(r)); namely, 



frequently observed 
quantity Q is characterized by a "local scaling" property 



((5p(r)) 2 8v{r)) = 2((p(x+r)) 2 V -v(x)). The above equa- 
tion is fundamental in relating the scaling of the density 
fluctuations in the inertial range to a flux- like quantity of 
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"density energy" analogous to the the case of a passive 
or active scalar. The single point value ((p(x)) 2 V • v(x)) 
should play the same role as "energy dissipation" does in 
the 4/5 law for the energy cascade in turbulence. This 
issue deserves more detailed study in the future. 

According to the multifractal theory of turbulence, the 
effect of saturation in the exponents z p and h p is equiv- 
alent to saying that the tail of the probability distribu- 
tion P[Sp(r)] should behave as ~ r z=o for any r. Thus 
we should expect that the functions r~ z °° P[Sp(r)] should 
collapse on the same "universal" distribution for r in the 
inertial range In figureQJwe show that this is exactly 
the case for ILES result. Let us note that, as shown in 
figurc[3l Zqo is larger than hoc ■ Using the multifractal the- 
ory, one can relate the saturation exponent to the fractal 
dimension of the "front-like" structures; i.e., D p = 3 — 
and Ds — 3 — hoo, where D p and D$ are the fractal di- 
mensions of the "front-like" structures for the density 
and entropy, respectively. Our findings show that Ds is 
larger than D p ; i.e., fronts in the entropy field are eas- 
ier to form than those in the density field. This result 
may not be surprising if we observe that large entropy 
fluctuations are produced by large pressure fluctuations 
or by large density fluctuations. Thus the fractal dimen- 
sion of the entropy fronts may be larger than those of 
density and pressure. The above argument implies that 
entropy is a more intermittent quantity than are density 
and pressure. The fact that the entropy satisfies a trans- 
port equation tells us it is a conserved quantity along La- 
grangian trajectories of the numerical simulations. This 
could suggest some connection between the existence of 
front-like structures and the behavior of inertial parti- 
cles in incompressible turbulence, where it is known that 
particles tend to form multifractal sets with correlation 
dimension as low as 2 [l8|, [3] . We feel that important re- 
sults can therefore be achieved by building a systematic 
"bridge of knowledge" between Lagrangian dynamics in 
compressible flows and entropy statistics in compressible 
flows. 

Let us summarize our results. Using a numerical 
simulation of inviscid homogeneous, isotropic weakly- 
compressible turbulence, we find that the scaling prop- 
erties of the velocity field in the inertial range are in 
excellent agreement with those observed in DNS of the 
Navier-Stokes equations This result supports the 

statement that the nature of the dissipation does not af- 
fect the statistical properties of the inertial range; i.e., 
turbulence is universal with respect to the dissipation 
mechanism. We confirm that transverse and longitu- 
dinal structure functions show different scaling proper- 
ties (up to this Reynolds number). We have also shown 
that, although almost no shock waves are produced in 
the simulation, the density fluctuations are characterized 
by "front-like" structures that determine the tail of the 
probability distribution of Sp(r). Accordingly, the scal- 
ing exponents of density and entropy structure functions, 



Dp and E p saturate at large p. 

The presence of front-like structures is well known 
in passive scalars advected by incompressible turbulent 
flows 1|| 17 1, which suggests that the density field be- 
haves as a passive-like scalar at small Mach numbers. 
The values of the saturation exponents for the density 
structure functions arc different from those observed in 
the 3D passive scalar case. This provides further confir- 
mation of the fact that the scaling properties of passive 
or "passive-like" quantities are not universal with respect 
to the properties of the advecting velocity field [2(j (here 
the correlation between density and the velocity is cer- 
tainly different from that observed for true passive scalar 
fields). The same argument may explain the difference 
between the scaling properties of the entropy we see in 
our simulation and those seen at larger Mach numbers in 

as. 
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